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Abstract We study the low temperature transport char- 
acteristics of a disordered metal in the presence of electron- 
electron interactions. We compare Hartree-Fock and dy- 
namical mean field theory (DMFT) calculations to in- 
vestigate the scattering processes of quasiparticles off 
the screened disorder potential and show that both the 
local and non-local (coming from long-ranged Friedel 
oscillations) contributions to the renormalized disorder 
potential are suppressed in strongly renormalized Fermi 
liquids. Our results provide one more example of the 
power of DMFT to include higher order terms left out 
by weak-coupling theories. 

Keywords Disorder • Strong correlations • Perturba- 
tion theory • Dynamical mean-field theory 
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1 Introduction 

Transport in dirty metals has been investigated for many 
years, with a substantial theoretical and experimen- 
tal understanding being achieved in the case of weak 
disorder and moderate electron-electron interaction [TJ 
|2]. However, much less is known about situations with 
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strong electronic correlations, where much of our under- 
standing relies on the application of numerical meth- 
ods like quantum Monte Carlo [3] or exact diagonal- 
ization [4], which are nevertheless severely limited in 
temperature range and system sizes. A more flexible 
approach to investigate the interplay of strong correla- 
tions and disorder is provided by the dynamical mean 
field theory (DMFT) [5 j. In its original formulation, 
the DMFT treatment of disordered systems does not 
include Anderson localization effects [6], a limitation 
which was quickly resolved by the introduction of the 
statistical DMFT (statDMFT) [7,8 . The statDMFT 
approach has already led to some novel effects like a 
strong disorder screening by interactions [6,9 , energy- 
resolved spatial inhomogeneities [9] , and the emergence 
of an Electronic Griffiths phase in the vicinity of the 
disordered Mott transition HOlfTTl. 



To partially elucidate the mechanism behind the 
rich physics uncovered through the statDMFT method, 
a recent work [12] provided analytical insights into the 
scattering off a weak random disorder potential in an 
otherwise uniform strongly interacting paramagnetic metal. 
While the analysis is most straightforward and trans- 
parent in this regime, this general issue is of key rele- 
vance also for the diffusive regime. Here, we revisit this 
problem, explicitly comparing the statDMFT results 
with those of the Hartree-Fock (HF) approximation. 
Our analytical results highlight the non-perturbative 
nature of the statDMFT and show that processes left 
out by HF generate vertex corrections to the impurity 
potential, which ultimately lead to enhanced screening 
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2 Model and methods 

We study the paramagnetic phase of the disordered 
Hubbard model 



icr (ij),<r 



(1) 



where Uj are the hopping matrix elements between near- 
est-neighbor sites, cj^c^) is the creation (annihilation) 
operator of an electron with spin projection a at site 
z, U is the on-site Hubbard repulsion, rii a = cl a Ci a is 
the number operator, and Si are the site energies (bare 
disorder potential). We consider here its paramagnetic 
phase at half- filling (chemical potential \i = U/2) and 
a particle-hole symmetric lattice. Below we discuss two 
different routes to treat this model. 



2.1 Hartree-Fock (HF) 

To solve the model ([!]), we first consider the HF ap- 
proach, as described, for example, in Refs. [13 , 14 . Here, 
we simply decouple the interaction term in as n^n^ p 
(n it ) nil + n n ( n ii) ~ ( n n) ( n ii)' We restrict ourselves 
to the paramagnetic solution, (n^) = (n^) = (n^), and 
the self-consistency condition is obtained from 



Ho 



(2) 



where T is the temperature, Gu (u n ) is the local part 
of the lattice Green's function, uo n are the Matsubara 
frequencies, Ho is the clean (ei = 0) and non-interacting 
(U = 0) lattice Hamiltonian, v is a site-diagonal matrix 
~ v i&iji whose entries 

■E i (P)-H = e i + U(n i )-lJL (3) 



are the renormalized HF disorder potential, and Ui (0) = 
Ei (uj) = U (rii) is the frequency-independent HF elec- 
tronic self-energy. We note that the HF approximation 
can be regarded as the static (weak-coupling) limit of 
the statDMFT and that it already contains one of its 
most important features: a self-energy which is local, 
but which varies from site to site reflecting spatial dis- 
order fluctuations. As we will show later, the statDMFT 
contains all the HF diagrams, but it also re-sums many 
higher order terms left out by HF. 

In general, we have to solve the self-consistency equa- 
tion in ([2| numerically. However, for a weak disorder 
potential (|^| <C D, where D is the bare half band 
width), we can expand it around the uniform solution, 
and, to leading order in the disorder potential, we have 

e(q) 



v(q) 



un n 



o 



(4) 





(a) 



(b) 



(c) 



Fig. 1 Diagrammatic representation of the electron- 
impurity interaction: (a) the bare electron-impurity inter- 
action vertex; (b) the first term of the Hartree-Fock the- 
ory, which is equivalent to the RPA resummation of "bub- 
bles"; (c) the first vertex correction to the electron-impurity 
interaction, absent in Hartree-Fock theory, but included 
in the statDMFT/slave boson solution. The inclusion of 
these and higher-order vertex corrections is essential for the 
phenomenon of perfect disorder screening. Here, the wavy 
line corresponds to the local on-site Hubbard-type electron- 
electron interaction. 



where e (q) is the inverse lattice Fourier transform of 
the disorder potential il q is the usual static Lind- 
hard polarization function [T5] . 

AsU —^0, the renormalized disorder potential reads 
Vi ~ ei + UIIijej, where 77^ is the lattice Fourier trans- 
form of il q , showing that the electrons scatter not only 
off the bare impurities, but also off the long-ranged po- 
tential generated by the Friedel oscillations (encoded 
in Ilij). This result contains another general feature of 
the statDMFT: the renormalized disorder potential ac- 
quires non-local terms, which are absent in the original 
DMFT formulation [6]. 

We can easily interpret Q in terms of the usual 
diagrammatic perturbation theory if we rewrite it as 

«(q) 1 



1 + U$ (q)7T q , 



e(q) krpa(<1) 
U 



u 



l-C/ilq «flPA(q)' 

where we have defined the dielectric function 



krpa(cl) = 1 -/777, 



(5) 
(6) 

(7) 



which in this approximation is given by the RPA ex- 
pression [T5] . It is then clear that the HF approxima- 
tion dresses the electron-impurity vertex by a "chain of 
bubbles" , as in the standard RPA screening theory [15] , 
Fig. 03). 

2.2 Slave-bosons (SB) 

To investigate the strongly correlated regime, we imple- 
ment the statDMFT using the slave-boson (SB) mean- 
field theory of Kotliar and Ruckenstein [16] (which is 
equivalent to the Gutzwiller variational approximation 



Mechanism for enhanced disordered screening in strongly correlated metals: local vs. nonlocal effects 



3 



[T7] ) as impurity solver jHHT]. This theory is mathe- 
matically equivalent to applying directly the original 
formulation of Kotliar and Ruckenstein to the Hubbard 
model in 0, as discussed in Ref. [9]. 

As in the previous HF calculation, we consider a 
weak disorder potential and expand the relevant mean- 
field equations [9j[TT] around their uniform solution. For 
particle-hole symmetry Zi = Zo + (sf) and we have, 
at T = [H] 



v(q) 



e(q) 



1 



un a 



o 



e(q)' 



where u = U/U c and 
1 



U 



1 



u. 



(8) 



(9) 



Eq. ([8| implies the following dielectric function 

^(q) = l-^ 2 -C/77 q . (10) 

The relation n (q) = 1/2 + 77 q v (q) also holds, as ex- 
pected. It is instructive to consider the limits of weak 
and strong interactions. If U <C J7 C , 

«Ss(q) ~ 1 - l/ilq = «HPA(q), (11) 

and we recover the HF result. In contrast, as U —> U c , 



Kss(q) -> 

leading to 



(i-e//e/ c ) c/- 1 [ii- 1 ],^,-. 



(12) 



(13) 



As the system approaches the Mott transition, the renor- 
malized disorder potential goes to zero at all lattice 
sites, a situation that was dubbed "perfect disorder 
screening" in Ref. [6 j. Its spatial structure is also very 
interesting, since V{ is just as non-local as for small U, 
but the non-local term is governed not by the Lindhard 
function, but by its inverse. The spatial structure of the 
charge disturbance in this limit is given by 



Srii 



[(1 - U/U c ) /U c ] 8 % + ((1 - U/U c f) . (14) 



Additional insight into these results can be obtained 
by noting that the second term on the right-hand side 



Thus, although the charge fluctuations are suppressed 
everywhere, its non-local part, coming from the Friedel 
oscillations, is much more strongly suppressed (O (1 — u) 3 ) 
and the electronic density is significantly disturbed only 
in the immediate vicinity of the impurities. The sup- 
pression of the slow spatial decay in Srti reflects the 
fundamental tendency of quasiparticles to become lo- 
calized as the system approaches the Mott insulator. 
Therefore, density fluctuations are healed very effec- 
tively in the strongly correlated limit. 



of Eq. ( 10 ) is unimportant both in the weakly and in 



the strongly correlated limits, cf. Eqs. (11) and (12). 



Neglecting this term in (10), we can follow the same 



procedure as in ([5| and rewrite (|8| as 



l + ^(q)i7< 



e(q) 



m (q) 



u 



i-un a 



(15) 



(16) 



The approach to Mott localization in this language can 
thus be described by the replacement U —> U, cf. Eqs. 
([5| and This replacement, in turn, may be viewed 
as a local field correction coming from vertex correc- 
tions in the polarization function [T5l [T8] (see Fig. ([1J) 
and the discussion in Section [3] below) . Close to the 
Mott transition U w U/(l — u) diverges and this strong 
correlation effect is seen to fall completely outside the 
scope of the HF theory. In fact, whereas HF predicts 
a gradual decrease of the dielectric function with in- 
creasing Z7, signaling the suppressed screening, see Eq. 
([7|, the statDMFT/SB approach predicts precisely the 
opposite: the dielectric function diverges as U —> U c , 



see Eq. (12), and the screening becomes asymptotically 
perfect. 



3 Beyond weak-coupling 

Comparing the renormalized disorder potential in Q 
and (J8|, we see that the interaction corrections left out 
by HF generate vertex corrections to the impurity po- 
tential, which are contained in the effective interaction 



U^(q) in (16). Since the HF approximation is the first 



term in expanding the electronic self-energy in U, we 
expand our solutions Q and ([8| in powers of Z7, in or- 
der to track down which terms are left out of HF. As 
we have seen, to first order in L 7 ", the statDMFT/SB 
and HF solutions agree. At second order, a difference 
emerges already 



(2) 
V HF 



(q) 



*(q) 



(2) 
V SB 



(q) 



*(q) 



u 2 nl 



u 2 u\ 



i 



(17) 
(18) 



To gain insight into the leading correction beyond 
HF, we combine the statDMFT procedure with usual 
perturbation theory. First, we recall that in the stat- 
DMFT approach the electronic self-energy is local, al- 
beit site-dependent. The only contribution to a local 
self-energy which is of second order in the interactions 
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is given by [19] 

4 2) (ico n ) = U 2 T J2 Gu } (i"n + iv n ) (iv n ) , (19) 



where 



n u (iv n ) = tJ2 4 0) K) of K - ^n) 



(20) 



is the local contribution of the dynamical Lindhard po- 
larization, calculated using the local Green's function 



4 0) (a*o 



■Wo)" 



with e» 



We are interested only in the leading uo n behavior of 



(19), since our statDMFT/SB approach is itself a low- 
energy one [20] . Ultimately, this low-energy approxima- 
tion provides a local Fermi-liquid description of the aux- 
iliary impurity problem [21] . For example, the uniform 

(2) 



contribution 
(l - 1/zj?) 



= 0) of 



lUJ n + U [jC- J, Willi I 1 

for a particle-hole symmetric lattice. 



(19) 
with ^ 



is 



Z { J («*>n) « 

1/Z< 2) ) ex -t/ 2 , 



We now expand (19) and (20) to linear order in the 



impurity potential. There are three identical contribu- 
tions, each with one of the three Green's function lines 
with an impurity vertex inserted in it, as shown in Fig 
(Jlj;). We focus on uf^ (0) since this defines the renor- 
malized disorder potential. Qualitatively, it is very easy 



to see how the extra terms in (18) are generated. Con- 



sider, for simplicity, that we estimate II a (iv n ) in ( |19| ) 
through the clean and static limit Ua (iv n ) « 77o (0) = 
— p (0) ~ —U~ x . In this case, we simply have u\ 2 ^ (0) ~ 
(U/U c ) 2 U c IIij£j, which has the same structure as the 



last term in (18). Based on these arguments, we stress 
that the difference which exists already at order U 2 



between (17) and (18) is an interaction-generated ver- 



tex correction of the electron-impurity vertex, which is 
absent in HF/RPA screening theory, but which is re- 
summed to all orders within statDMFT/SB. 



4 Conclusions 

We presented a detailed analytical calculation of the 
effects of weak disorder scattering in a correlated host. 
Comparing the results obtained within HF and stat- 
DMFT, we highlighted the fact that statDMFT incor- 
porates important vertex corrections to all orders, a 
task which is difficult, or more likely, even impossible to 
perform using weak-coupling diagrammatic approaches. 
A physical consequence of the inclusion of these vertex 
corrections is the phenomenon of disorder screening by 
interactions. 

An analogous example of this dichotomy can be ob- 
served in the familiar Migdal-Eliashberg (ME) strong 



coupling theory describing the electron-phonon prob- 
lem [15]. Indeed, the ME theory neglects the momen- 
tum dependence of the electronic self-energy and may 
thus be regarded as a weak-coupling approximation to 
DMFT. Like the HF/RPA screening described above, 
it also neglects vertex corrections. The full DMFT solu- 
tion, however, not only contains all the ME diagrams, 
but it also re-sums many higher order terms left out by 
the ME approach, including vertex corrections [22], in 
close analogy with the statDMFT treatment of disorder 
and interactions. In both models, these strong coupling 
effects reflect non-perturbative Kondo-like processes [6j 
[22] . which cannot be described by weak-coupling ap- 
proaches. 
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